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ABSTRACT 

We present a counter example which shows the violation of the S-matrix 
factorization in the massive Thirring model. This is done by solving the PBC 
equations of the massive Thirring model exactly but numerically. The violation 
of the S-matrix factorization is related to the fact that the crossing symmetry 
and the factorization do not commute with each other. This confirms that the 
soliton antisoliton S-matrix factorization picture of the sine-Gordon model is 
semiclassical and does not lead to a full quantization procedure of the massive 
Thirring model. 
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1. Introduction 

Integrable field theories in 1 + 1 dimensions have been studied in various 
fields of research. The Heisenberg xxx model is solved by the Bethe ansatz 
technique [1]. The solution is confirmed most elegantly by the Yang-Baxter 
equation [2-4]. Also, the Heisenberg xxx model with spin | is equivalent to 
the Thirring model which is solved exactly in the continuum field theory. 
Since the Yang-Baxter equation corresponds to the S-matrix factorization, 
there is, therefore, no question about the S-matrix factorization in this case. 

Also, the Heisenberg xyz model with spin | has been studied quite ex- 
tensively [5-9] . Since this model is proved to be equivalent to the massive 
Thirring model in the continuum limit [10], it should be interesting to study 
these models from different points of view. 

The S-matrix factorization is also assumed for the massive Thirring model 
or sine-Gordon field theory [11]. With this ansatz of S-matrix factorization, 
the spectrum of the sine-Gordon model has been obtained and is found to 
be consistent with the spectrum obtained by the WKB method [12-14]. 
Therefore, it has been beheved that the factorization of the S-matrix for 
the massive Thirring model holds exactly since the bound state spectrum 
obtained by the WKB method was considered to be exact. However, the 
recent investigations of the bound state problem of the massive Thirring 
model from the light cone method as well as from the Bethe ansatz tech- 
nique have shown that there is only one bound state, on the contrary to the 
WKB result [15-18]. Indeed, it is proved that the WKB result of the bound 
state spectrum for the massive Thirring model is not exact [15-17]. This 
suggests that we should reexamine the S-matrix factorization for the mas- 
sive Thirring model since the S-matrix factorization has also been believed 
to hold in an exact fashion. 

Until now, this factorization ansatz of the S-matrix for the massive Thirring 
model has never been questioned seriously. In fact, the factorization of the 
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S-matrix for the particle particle scattering holds exactly. This can be easily 
seen by using the Bethe ansatz wave functions. However, it is nontrivial 
to obtain the S-matrix factorization for the particle-antiparticle scattering 
case, as we will see it later. 

Now, the point is that the factorization of the S-matrix has been obtained 
for the sine-Gordon field theory with sohton anti-sohton pictures. There- 
fore, in this case, one docs not have to worry about the antiparticle property. 
However, the soliton and antisoliton are only semiclassical objects and there 
is no clear way to quantize it. This is related to the fact that the soliton is 
a solution of the field equation, but not the field itself. 

Therefore, even though one can easily write down the particle antiparticle 
scattering S-matrix with the factorization ansatz starting from the particle 
particle S-matrix factorization property, the identification of the soliton 
and antisoliton as the fermion and antifermion cannot be justified quantum 
mechanically. 

In this paper, we present a counter example which shows the violation of 
the S-matrix factorization for particle hole scattering case in the massive 
Thirring model. This is based on the Bethe ansatz solutions which we 
obtain exactly by numerically solving the PBC equations. The structure of 
the Bethe ansatz solutions is quite similar to the free fermion field theory. 
The ultraviolet cutoff A is determined by the fermion number when we 
put the theory into the box with the length of L. In free field theory, the 
A is given as 



while, in the massive Thirring model, A can be obtained as the function of 
the N/ L after we solve the PBC equations properly. Any physical observ- 
ables can be obtained by letting the N and L infinity. Here, we note that 
this is the simplest but the best way to define the field theory. 

In order to discuss the S-matrix factorization, we have to determine the 
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vacuum and then create particle hole states as long as we are interested 
in the charge zero sector. Since the number operator commutes with the 
hamiltonian of the massive Thirring model, we are only concerned with the 
charge zero sector. 

Here, we show that the S-matrix factorization does not commute with the 
crossing symmetry. Since the lagrangian of the massive Thirring model 
is invariant under the charge conjugation, one tends to believe that the 
crossing symmetry should always hold. Indeed, the crossing symmetry it- 
self holds. However, one must be careful for the order of operations in 
quantum mechanics. There is no guarantee that one can make use of the 
crossing symmetry for the S-matrix of the particle hole scattering with the 
factorization properties together. 

In order to check the validity of the S-matrix factorization, we have to solve 
the equations of the periodic boundary conditions (PBC) and construct the 
vacuum. However, this is not sufficient if we want to deal with the particle 
hole scattering. We have to solve the PBC equations for n— particle n— hole 
states. 

The important point is that the PBC equations for n— particle n— hole states 
are different not only from those for the vacuum states but also from each 
other. Therefore, we have to first solve these PBC equations to determine 
the rapidities of the vacuum as well as those of the n— particle n— hole 
states. 

It is always an interesting question how one can construct the particle hole 
state under the situation where the rapidity distribution of the particle hole 
states are different from those of the vacuum state. In other words, how 
can we find or identify the rapidity of the hole state if the rapidities of 
the negative energy particles for the n— particle n— hole states are different 
from those of the vacuum state ? Later, we will see that we can find the 
rapidity of the hole state in the limit of infinitely large N (the number of 
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particle) and L (the box size). In this hmit, we can maintain the particle 
hole picture. 

Here, a question arises. Can we also prove the factorization of the S-matrix 
even for the particle hole scattering case ? As can be seen from the above 
statement, the hole states carry the information of all the negative energy 
particles since the rapidities of the particle hole state are different from the 
vacuum. Therefore, it is highly nontrivial that the S-matrix for the particle 
hole scattering can be factorized. 

Indeed, we will show below that the S-matrix for the particle hole scattering 
is not factorizable since the factorization and the limit of N and L infinity 
do not commute with each other. We stress here that, if we show any 
examples of the violation of the S-matrix factorization, this is sufficient 
that one cannot rely on the factorization ansatz. On the contrary, if we 
wanted to show the validity of the S-matrix factorization, then it would 
have been a very hard work. 

The violation of the factorization ansatz indicates that the bound state 
spectrum constructed from the S-matrix factorization cannot be justified 
any more. As mentioned above, this is consistent with recent calculations 
on the massive Thirring model, which show that the Bethe ansatz PBC 
equations give only one bound state, on the contrary to the semiclassical 
result and that of the Bethe ansatz technique with the string hypotheses. In 
fact, it is shown that there is no string-like solution that satisfies the PBC 
equations if they are solved exactly for the particle hole configurations [15]. 
Also, the light cone prescription of the massive Thirring model shows that 
there is only one bound state [16,17]. This indicates that the S-matrix 
factorization must correspond to the semiclassical approximation. This 
point is indeed shown in this paper, that is, the S-matrix factorization can 
be justified only when one neglects the commutability of the operators. This 
means that the sohton antisoliton picture is indeed semiclassical. Therefore, 
in order to make a correct correspondence between the soliton and fermion 
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in a fully quantized sense, we should not rely on the S-matrix factorizations. 

This paper is organized as follows. In the next section, we will briefly de- 
scribe Zamolodchikov's S-matrix factorization method. Then, in section 3, 
we discuss the Bethe ansatz solution for the massive Thirring model. In 
section 4, we present the numerical solutions of the Bethe ansatz PBC equa- 
tions in order to define the vacuum as well as the n— partcile n— hole states. 
Then, section 5 will treat the S-matrix for the particle hole scattering. In 
particular, we will present the numerical evidences that the S-matrix of the 
particle hole scattering is not factorizable. Section 6 will summarize what 
we have understood from this work. 



2. Zamolodchikov's S-matrix factorization 

In integrable field theories, there are infinite conservation laws which in- 
duce strong selection rules [11]. For example, there is no particle creation 
or annihilation in the scattering process. Also, there is no change of the 
momentum before and after the scattering process. This means that the 
scattering is always elastic. 

The S-matrix for the two body scattering case can be defined as 

S\Ai{ai)Aj{aj) >i„= ^ S^j{ai, aj)\Ak{ai)Ai{aj) >out (2.1) 

k,l 

where Ai{a)^s describe the particle states with the rapidity a. 
In the same way, we can consider the N body S-matrix 

h—3N 

(2.2) 

The basic assumption of the S-matrix factorization is that this S-matrix 
can be written by the product of the two body S-matrices. 
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For simplicity, we consider the three body case. One assumes that the 
S-matrix of the three body scattering can be written as 

SlT- E ^^'^^(«i,«2)5^,\"(a3,«i)5S:,3(«2,«3) (2.3) 

This assumption can be checked by making use of the Bethe ansatz solutions 
to the massive Thirring model for the particle-particle scattering case. 

In this three body scattering case, we can consider the following two differ- 
ent processes. In the first case, the particle 1 scatters with the particle 2, 
and then the particle 1 scatters with the particle 3, and finally the particle 
3 scatters with the particle 2. On the other hand, in the second case, the 
particle 2 first scatters with the particle 3, and then the particle 3 scatters 
with the particle 1, and finally the particle 1 scatters with the particle 2. 

Since the final state is the same between the above two processes, the two 
scattering events should give the same scattering process. Therefore, we 
obtain 

Pl,P2,P3 

= E S]l''^^{a,,as)Slll{a,,a^)S^^^^{a,,a,) (2.4) 

91,92,93 

This is the Yang-Baxter equation. 

Now, Zamolodchikov further assumes that the S-matrix should satisfy the 
unitarity, the analyticity and the crossing symmetry. The unitarity reads 

h ,32 

The analyticity is written as 

S\a)=S{-a*). (2.6) 

Finally, the crossing symmetry can be stated as 

SS{a)^S^{t7r-a) (2.7) 

where i and I denote the anti-particles of i and / states. Here, the an- 
tisoliton of the sine-Gordon model is identified as the antifermion of the 
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massive Thirring model. Together with the above assumptions, Zamolod- 
chikov obtained the functional equation with the particles and antiparticles 
included. 

However, for the massive Thirring model, it is nontrivial to define the an- 
tiparticle state in the nonperturbative sense. This is because one has to 
first determine the vacuum and then construct the particle hole states. 
The problem here is that the rapidity distributions for the vacuum and the 
particle hole states are different from each other since the periodic bound- 
ary condition equations are different. Therefore, it is nontrivial to make 
hole states referring to the vacuum rapidity distribution. 

The question we want to address here is whether this crossing symmetry 
(eq. (2.7)) can be satisfied together with the factorization of the S-matrix 
or not. 

It is obvious that the crossing symmetry itself should hold since the La- 
grangian of the massive Thirring model possesses this symmetry. Indeed, 
this can also be proved by the Bethe ansatz solutions. 

However, it is highly nontrivial whether the factorization for the particle 
hole scattering S-matrix can also hold or not. This is what we want to check 
in this paper and we prove numerically that the particle hole S-matrices do 
not satisfy the factorization. 



3. Bethe ansatz solution of massive Thirring model 

The massive Thirring model is a 1-1-1 dimensional field theory with current 
current interactions [19]. Its lagrangian density can be written as 

L = i,{i^^d^ - mo)V' - \9^f3^ (3.1) 
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where the fermion current is written as 

jn -■ V'TmV' : ■ (3.2) 
Choosing a basis where 75 is diagonal, the hamiltonian is written as 
H = J dx 



d d 



(3.3) 

Now, we define the number operator N as 



N 



j dxil}^. (3.4) 



This number operator commutes with H. Therefore, when we construct 
physical states, we must always consider physical quantities with the same 
particle number N as the vacuum. For different particle number state, the 
vacuum is different and thus the model space itself is different. 

(a) Bethe ansatz wave functions 

The hamiltonian eq.(3.3) can be diagonalized by the Bethe ansatz 
wave functions. The Bethe ansatz wave function ^(xi, xjv) for 
particles can be written as [13-14] 

^(xi, ...,a;Ar) = exp(imo^a;jSinh/3i) J| [I + i\{(5i, (5j)€{xi - Xj)] 

l<i<j<N 

(3.5) 

where /3i is related to the momentum ki and the energy Ei of the i-th 
particle as 

ki — mo sinh (3.6a) 

Ei = mo cosh Pi. (3.66) 

Here, /3i's are complex variables. 

e{x) is a step function and is defined as 

r -1 x<o 



e{x) — < 



(3.7) 

1 x>0. 



X{Pi,Pj) is related to the phase shift function (f){Pi — Pj) as 

expMA-/3,)) = ^— 
The phase shift function — pj) can be exphcitly written as 



(j){Pi - Pj) = -2 tan"^ 



(3.8) 



(3.9) 



In this case, the eigenvalue equation becomes 



TV 



(3.10) 



H I Pi...Pn >= C^mocosh Pi) \ Pi...Pn > 

i=l 

where | Pi-.-Pn > is related to ^{xi, ...,xn) as 

N 

I Pi...Pn >= / dxi...dxiv*(xi, xn) n i^H^i^ A) I > . (3.11) 

i=l 

Also, ip{x,P) can be written in terms of ipi{x) and ■02(2^) as, 

V^(x,/3) = e2V'i(a;) + e"2V'2(a:^)- (3.12) 

Prom the definition of the rapidity variable /9j's, one sees that for posi- 
tive energy particles, PiS are real while for negative energy particles. Pi 
takes the form iir — ai where cti's are real. Therefore, in what follows, 
we denote the positive energy particle rapidity by Pi and the negative 
energy particle rapidity by ctj. 



(b) Periodic Boundary Conditions 

The Bethe ansatz wave functions satisfy the eigenvalue equation [eq.(3.10)]. 
However, they still do not have proper boundary conditions. The sim- 
plest way to define field theoretical models is to put the theory in a 
box of length L and impose periodic boundary conditions (PBC) on 
the states. 

Therefore, we demand that "^{xi, ..,xn) be periodic in each argument 
Xi. This gives the boundary condition 

^{Xi = 0) = ^{Xi = L). (3.13) 
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This leads to the following PBC equations, 

exp{imoLsmhPi) = exp(-i^0(/3j - Pj)). (3.14) 

i 

Taking the logarithm of eq.(3.14), we obtain 

moL sinh /3j = 2Tmi — ^ 4>{(3i — j3j) (3.15) 

j 

where n^'s are integers. These are equations which we should now 
solve. 



4. Numerical Solutions 

The parameters we have in the PBC equations are the box length L and the 
particle number N. In field theory, we introduce the ultraviolet cutoff pa- 
rameter. In the PBC equations, the cutoff parameter is the particle number 
N. This is quite similar to the free field theory where the cutoff parame- 
ter is defined by the particle number N. Once we determine the particle 
number N and the box length L, then we determine all the rapidity values 
necessary to obtain any physical quantity. In this sense, the Bethe ansatz 
solution can be well defined as a field theory. Any physical observables can 
be obtained by letting the N and L infinity. 

There is also an important quantity which is the density of the system p. 
It becomes 

(4.1a) 

Here, the system is fully characterized by the density p. For later conve- 
nience, we define the effective density po as 
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where Lq and A^o are defined as Lq = moL and A^o — |(-^~ 1); respectively. 

It is important to note that the beta function of the massive Thirring model 
vanishes to all orders [20], and thus there is no need to worry about the 
renormalization of the coupling constant. However, one has to be careful 
for the coupling constant normalization arising from the way one makes 
the current regularization. Here, we employ the Schwinger's normalization 
throughout this paper. 

In what follows, we solve the PBC equations numerically in the same man- 
ner as in ref.[15]. The numerical method to solve the PBC equations is 
explained in detail in ref.[15]. We note that the errors involved in the nu- 
merical calculations are very small. For example, the numerical uncertainty 
for the rapidity values will appear at the 10~^ level, and therefore, this does 
not generate any problem in the present discussion. 



(a) Vacuum state 

The PBC equations for the vacuum which is filled with negative energy 
particles { Pi — iir — ai ), become 



. 27rni 2 ^ _^ 
smh ctj = — 1~ 

L/n Lin ■_/• 



1 , 1 

-c/0 tanh - 



Qq tanh -(ctj — aj) 



(z = l,..,7V) 
(4.2) 

where rij runs as 

= 0,±l,±2,...,±iVo. 

There is no ambiguity to determine the vacuum. The vacuum is indeed 
constructed uniquely [15]. 

(b) Ip — Ih configuration 

One particle-one hole (Ip — Ih) state can be made by taking out one 
negative energy particle (the riig particle) and putting it into a positive 
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energy state. In this case, the PBC equations become 



Hi ^ n. 



sinh ai — 



— —tan 



2 V- -1 [1 ,1/ 
-— 2^ tan -goteinh-{ai- aj) 

Lin I - ^ 



(4.3a) 



rii 



rii, 



sinh p. 



^ + T- E tan-^ \lgo coth l(Ao + #)^b) 



where can be a complex variable as long as it can satisfy eqs.(4.3). 

It is important to note that the momentum allowed for the positive 
energy state must be determined by the PBC equations. Also, the 
momenta occupied by the negative energy particles are different from 
the vacuum case as long as we keep the values of N and L finite, as 
can be seen from eqs.(4.2) and (4.3). 

Note that the lowest configuration one can consider is the case in which 
one takes out = particle and puts it into the positive energy 
state. This must be the first excited state since it has a symmetry of 
ctj = — Q;_j. Indeed, as discussed in ref.[15], this state corresponds to 
the only bound state (the boson state) in this model. 
Next, we consider the following configurations in which we take out 
rii — il) i2, .. particles and put them into the positive energy state. 
These configurations of one particle-one hole states turn out to be the 
scattering states [15]. 



(c) 2p — 2h configurations 

In the same way as above, we can make two particle-two hole (2p — 2h) 
states. Here, we take out the and rijj particles and put them into 
positive energy states. The PBC equations for the two particle-two 
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hole states become 



2Tmi 2 

smn ctj = — — tan 

Lq Lq 



^5'ocoth^(aj + /3jJ 



2 



tan" 



-1 



^90 coth^(ai + 



^ tan ^ ^5fotanh^(Q;i - aj) 

l^n ..... \_Z Z 



(4.4a) 



Ui = n,; 



sinh /3, 



27mi, 2 1 
— + — tan-i 



^^otanh^(/3i, - pi. 



2 rl 1 

+ — tan"^ -5roCoth-(/3j^ +0;^) 



(4.4b) 



n,- = n, 



«2 



sinh — 



2Tmi^ 2 _i 
+ — tan 



^t/otanh^(/3i2 - /3j J 



2 rl 1 



(4.4c) 



All the configurations one can construct for the two particle two hole 
states turn out to be the scattering states in this model [15]. 



5. Factorization of S-matrix 

Since the vacuum state is denoted by the rapidities ai, .., ctjv, 

|vac >= ctAT > 

we can define the S-matrix of the vacuum to vacuum transition. In this 
case, we can write it as 

N 

< vac|5'|vac >= JJ SQ{ai, aj) (5.1) 
14 



where <S'o(Q;i, aj) denotes the two body S-matrix between i and j particles 
to make a transition from Xi < Xj to Xi > Xj states. It can be written as 

So{ai, aj) = exp {i(f){ai - aj)) 

where (f){ai — aj) can be written exphcitly as 

(5.2) 



(j){ai — aj) = 2 tan ^ -qq tanh -(ojj — aj) 



.2" 2 

Note that these particles here are negative energy ones. Obviously, the 
vacuum to vacuum transition must be unity, 

< vac|5'|vac >= 1. (5.3) 

Therefore, we have always the constraint 

N 

l[So{ai,aj)^l. (5.4) 



We also give a correspondence between 5'o(ai, aj) as defined here and SP^{ai, aj) 
as defined in section 2. 



Sppi^i^ «j) = ^o{(^i, c^j) (5.5) 

where pp denotes the particle particle scattering. If it is for particle hole 
scattering case, the S-matrix can be written as SpJ^ai, Pj) as will be soon 
explained below. It should be noted that, in the massive Thirring model, 
the states are completely specified by the rapidity variables a with the index 
of particle or hole. These are the quantum numbers that can specify the 
states. 



(a) S-matrix 

Now, we first define the one particle one hole state. This can be 
denoted as 

\lplh >= |/5i,4,-,«]v >= > (5.6) 
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where is the rapidity of the hole state when N and L become 
infinity. But for the moment, it is only symbolically written since aj's 
differ from a^'s of the vacuum state. 

Therefore, the S-matrix for the Ip — Ih state can be written as 

TV TV 

'5ipi/i(/3i,a2,--,«jv) = n Soioii,al)Y[So{Pi,al). (5.7) 

i>j>l k=2 

Here, we note that, in the large N and L limits, we find that 
as will be shown later numerically in Table 1. 

In what follows, we only discuss those particle hole scattering processes 
where the particle states are always involved. We do not consider here 
particle particle scattering or hole hole scattering processes. 
Wc note that this particular example is sufficient since we are only 
interested in finding an example of the violation of the factorization. 

Therefore, we define the S-matrix for the particle-hole scattering in 
the following way, 

5;^(A,/3f) = nWi,4)- (5-8) 

k=2 

Next, we want to define the two particle two hole state. This can be 
denoted as 

\2p2h >= P2, al .., >= /3f , P2, > ■ (5.9) 

Again, and P2 written symbolically. 

Therefore, the S-matrix for the 2p — 2p state can be written as 

S2p2h{Pl, P2, 0(1, .., a]^) 
JV TV TV 

^So{Pi,P2) n So{ala^^)l[So{Pi,al)l[So{P2,aj). (5.10) 

i>j>2 k=3 1=3 
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Here again, we note that, in the large N and L hmits, we find that 

In the same way as SpJ^, we define the 2p — 2h scattering process where 
the particle states are involved. 

N N 

SX'M, P2, P^) ^ n SoiPi, 4) n So{P2, aj). (5.11) 

fe=3 1=3 

Now that the aj and q;| approach to the ctj for the large N and L 
limits, one may think that one can replace the aj and aj by the ai. We 
denote these S-matrices by /3?){v) and SXhiPi, Pi, P2, P^){v). 

Therefore, if we take the large N and L limits at this stage without 
care for the difference between ai, a} and aj, then we obtain that the 
S-matrix Sp^pl{(3i, P2, P2)(y) for the 2p — 2h state can be factorized 
into Ip — Ih S-matrices. That is. 




(5.12) 



Eq.(5.12) shows that the S-matrix of the two— particle two— hole scat- 
tering seems indeed factorized. This is the main reason why people 
believe that the factorization ansatz for the particle hole scattering 
should hold in the same way as the massless case or in other words as 
the particle particle scattering case. 

Based on the factorization ansatz of the S-matrix, Zamolodchikov and 
Zamolodchikov obtained the Yang-Baxter type functional equation for 
the S-matrix which can determine the shape of the S-matrix itself 
by the simple algebraic manner [11]. However, the derivation of the 
functional equation implicitly assumes the factorization of the particle 
hole scattering S-matrix. Even though the factorization of the particle 
hole S-matrix is valid in terms of their rapidity variables, it is nontrivial 
to show that this factorization can hold even for the common rapidity 
variables with the vacuum. 
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Therefore, one must be careful for the treatment of the rapidity vari- 
ables when one calculates the S-matrix. In particular, one should 
carefully treat the order of the operations when one makes the large 
N and L limit in quantum mechanics. The large limits of N and L 
should be taken as late as possible. 



(b) Large L and N limits 

In what follows, we show that, in the large N and L limits, eq.(5.12) 
does not hold. In order to see it, we first rewrite the S-matrix of the 
Ip — Ih state in terms of the vacuum rapidities. 

TV 



SpHWuP'^) = n SoiPi^ «fe + Sk) (5.13) 

k=2 



where 5i is defined as the difference between ctj and aj, that is. Si 



Since Si is quite small for the large values of N and L, we can expand 
So{(3i, aj + Sj) in terms of Si. Here, we expand the phase shift function 
+ aj + Sj — in) 



So{Pi,aj + Sj) = exp 



d 

i(j){Pi + aj — in)) + iSj-^4){l3i aj — in) 



(5.14) 



The first terms of eqs.(5.14) correspond to the S-matrices which are 
factorizable. 

Now, the important point is that the Si itself becomes zero as the values 
of and L become infinity but the summation of Si stays finite in this 
limit. We will discuss it in detail below. Therefore, it is important to 
keep the second term for the S-matrix evaluation. 

In the same way, we rewrite the S-matrix of the 2p — 2h state. 

'5p%t(/3i, Pi, P2, P^) = n SoiPi, + ek) n SoiP2, ai + q) (5.15) 

fe=3 1=3 
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where ej is defined as ej = q;| — a j. Since is very small for the large 
values of N and L, we can expand 5'o(/5i, + e^) and 80(^2, Oii + ei) 
in terms of in the same way as eqs.(5.14) 



So{Pi,ak + ek) =exp 



d 

i(f){Pi + ak- in) + iefe- — + at - in) 

dak 



i (9^ 



(5.16a) 



So{P2,Oii + e;) = exp 



d 



i(j){P2 + in) + iei——(j){P2 + ai - in) 

oai 



+ 2^/^^<^(/^2 + ai- in) + ... 



(5.166) 



We can also check later that the itself becomes zero when and L 
become infinity. However, the summation of is finite and therefore, 
we have to keep them when we want to show the S-matrix factorization. 

Now, we rewrite eqs.(5.13) and (5.15) in terms of eqs.(5.14) and (5.16). 
As mentioned above, we are only concerned with the particle hole 
scattering S matrix factorization. 

Therefore, eqs.(5.13) and (5.15) can be written as 

TV 



'5pt(/3i,/3i')=exp 



J2 HPi + ttm - in) 



m=2 



X exp 



TV 



k=2 \ 



d(t){(3i + ak-in) 1 ^ d'^HPi + "fc - in) 



dal 



(5.17) 



sXh{Pi,Pi,P2,f3^)^ew 



TV 



TV 



(f){pi + am-in) + iJ2 </'(/^2 + - ^tt) 



m=3 



m=3 



X exp 



N 



. ^ , d(j){f3i + ak-in) 1 ^ d'^cjy^Pi + au - in) 

^ + o^fc- 



fc=3 



dak ' dal 



/ (9(/)(/32 + ak-in) 1 2 d'^(t){P2 + ak-in^ 

+ ^ 2^ efc ^ n^k I + ••• 



k=3 



dai 



dal 



. (5.18) 
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Clearly from eqs.(5.17) and (5.18), the S-matrix factorization cannot 
hold unless the exponential terms except the first one vanish all to- 
gether. 

The important point is that J2 contains all the information of the 
Ip—lh state and therefore depends on (3i and which are the quan- 
tum numbers that characterize the Ip — lh state. In the same manner, 
the summation J2^k depends on Pi, (3^, P2 and which characterize 
the 2p — 2h state. 

(c) Numerical results 

Now, we want to present the numerical check of the various quantities 
which appear in eqs. (5.17) and (5.18). Here, we first solve numerically 
the PBC equations (4.2), (4.3) and (4.4). We vary the number of 
particle and the box size Lq. 

First, we want to check how the rapidity values of the Ip—lh and 
2p — 2h configurations converge into those of the vacuum state when 
N becomes large. In Table 1, we plot several values of the rapidity as 
the function of N with the fixed value of the density po- This clearly 
shows that the rapidities of the Ip — lh and 2p — 2h states approach 
to those of the vacuum state when N becomes very large. Therefore, 
eq.(5.12) indeed holds. However, we must be careful when we treat any 
physical quantities which depend on the sum of the rapidity difference 
between the vacuum state and the Ip — Ih or 2p — 2h configurations. 
The important point is that, even though each rapidity of the Ip — Ih 
and 2p — 2h states approaches to that of the vacuum, the sum of the 
rapidity differences stays finite. 

In order to see the effect mentioned above, we define the following 
quantities 



N 



d(f){Pi + Qfc - iir) 
dak 



(5.19a) 



k=2 
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E,Mn.) - E^MPu P',) = I E 5,2^W±|^ (5.196) 

^ k=2 ^^k 

k=3 "^k 

(5.19c) 

F (r, r,\--F (R R Rh^_^^ . d^(l){(3i + ak - in) 

i^2p2h\P'l,n2) = ^2p2hyPl, P2,Pl,P2) — 2 2^ ^^2 

fe=3 k 

(5.19rf) 

where j3i and /32 correspond to the rapidity values for the ni and n2 
states. 

3i+Qfc-»^) oTi,^ 9^</'(/3i+«fc-i7r) 

eq.(3.9) as 

+ Qijt - i7r) go 



Here ""^^^X?' ^nd '^^^1^+^'^-^^^ can be written exphcitly using 



(5.20) 



dau 2 | + (i + ^)sinh^i(^^ + Q;,) 

2 

a20(/3i + Qfc - i7r) _ go (1 + ^) sinh \{l3i + ctfc) cosh |(/ji + ak) 
" 2 (| + (i + |)3inh2i(/?, + a,))' 

(5.21) 

Since we have solved the PBC equations numerically, we know all the 
rapidity values, and thus we can calculate eqs.(5.19) together with 
eq.(5.20) and eq.(5.21). 

In Table 2, we plot the values of Dipi/i, £^ipi/i, -D2p2/i and £^2p2fe as the 
function of N and Lq for the lowest scattering state. As can be seen, 
the calculated values of Dipih and D2p2;i do not depend very much on 
the N and Lq as long as we keep the same density po- 
However, as can be seen from Table 2, the Ei^ih and the £^2p2/i decrease 
as N increases. They behave as 

Eipih{ni) ~ ^ 
E2p2h{ni,n2) ~ 

Therefore, these quantities do not survive at the large N limit. 
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This clearly shows that the sum of the rapidities times the derivative 
of the S-matrix is finite even though each rapidity value becomes zero 
when and L become infinity. Therefore, we prove that these second 
terms do not vanish. 

In order to see more clearly the above effects, we rewrite the Sp^{(3i, (5^) 
and Splp^{f3i,/3^,f32,l3^) in terms of Dipih{ni) and -D2p2?»(ni, 712). 

SplWi, = S^liP,, (3^\v) exp [iD.Mni)] (5.22) 

'S'pfepft(/5l,/5{',/52,/52) = S^l{Pl,Pl)(V)Sphi.P2,P2){V) 

X exp [iD2p2h{n\, 712) + iD2p2h{n2, rii)] (5.23) 

Now, it is obvious that the -Dipife(ni) and D2p2h{ni, 77,2) arc quite differ- 
ent from each other since they carry the information of the many body 
nature of the particle hole states. Indeed, the D2p2h{ni,n2) depends 
on /32, Pi and Therefore, the factorization of the particle hole 
S-matrix cannot be justified. In fact, as we see below, the numerical 
results confirm the difference between the Dipih{ni) and -D2p2/i(^i, ^2) 



(d) Field theory limit 

Now, we should take the field theory limit since this is the only way 
to compare our results with Zamolodchikov's factorization ansatz. For 
the field theory limit, we have to let p — > 00 [15]. In the case of the 
bound state problem, when we let p 00, then we should take the 
limit of mo 0, keeping the excited state energy finite. In order to 
see it more concretely, we repeat here how it is done for the bound 
state problem. The excited state energy AE is written as [15] 

AE = ra„ {a {£-)') 

where A and B are some constants which depend on the coupling 
constant. Here, a is a constant with positive value which is smaller 
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than unity. In this case, when we make p — > oo, keeping /S.E finite, 
we can make a fine tuning of mo such that 

mj-"p" = finite. 

This means that we should let mo — > 0. 

However, in the present case, the situation is a bit different and even 
simpler since the correction factor itself does not have any dimensions. 
Therefore, the field theory limit is just that we let the po infinity. In 
fig.l, we show the calculated absolute values of Dipih and D2p2h as the 
function of Pq^. In this case, they can be well parameterized by the 
following functions as 

DiM^i) = C(^Vo"' + D^"-^ (5.24a) 

^2p2/.(ni, ris) = C(') exp + L'^^) (5.246) 

where C^, D^^ and H2 are some constants which depend on the cou- 
pling constant and the rapidity values of the particle hole states. 

As can be seen, when the po goes to infinity, then the Dipih and D2p2h 
approach to some finite numbers. Namely, they become 

DipihM ^ D^^^ (5-25a) 
D2p2h{nu 712) ^ C(') + (5.256) 

This means that we have made the field theory limit in a correct way. 
The values of C('\ D^'^^ and K2 for several cases of the particle hole 
excited states are shown in Table 3. 

From the fig.l, we know that the breaking of the factorization of the 
S- matrix is the order of — D2p2h\ which is smaller than unity. 

We do not know whether this value of — D2p2h\ can be said to 

be large or small. In any case, the S-matrix factorization is violated 
at the elementary level. 
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Here, we again note that the S-matrix factorization for the particle 
particle scattering must hold. This is essentially the same as the vac- 
uum case. It is interesting to observe that the energy of the vacuum 
state can be obtained analytically, though the vacuum energy alone is 
not physically very meaningful. 



6. Conclusions 

We have presented a numerical proof that the S-matrix factorization for 
the particle hole scattering in the massive Thirring model is not satisfied. 
This is mainly because the factorization of the S-matrix and the crossing 
symmetry do not commute with each other. 

It is always simpler to find a counter example of the violation of the S- 
matrix factorization than to prove the validity of the factorization. This 
is clear since we only have to find out any type of example which shows 
the violation of the factorization. This is indeed the point we have worked 
out in this paper. The present calculations present a convincing evidence 
that the S-matrix factorization for the particle hole scattering does not hold 
exactly. 

Here, we want to examine the implication and consequence of this result. 
As stressed in this paper, it is important to realize that the S-matrix factor- 
ization should hold for the particle particle scattering even for the massive 
Thirring model. This can be easily seen from the Bethe ansatz solution 
which is eq.(3.5). Indeed, if one looks at the vacuum to vacuum transition 
eq.(5.1), then one sees that the S-matrix is factorized into the product of 
the two body S-matrices. This is the consequence of the particle particle 
scattering. Clearly, the massless case must have this nice feature of the 
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S-matrix factorization since, in this case, one can treat all the scattering 
processes as two species of fermion scattering without going to the particle 
hole scattering. Therefore, those models which are equivalent to the six 
vertex models should receive no effects from the present investigation. 

However, those models which are equivalent to the eight vertex models 
should be carefully treated if one is concerned with the bound state spec- 
trum with S-matrix factorization ansatz. At least, we now know that the 
massive Thirring model docs not have the S-matrix factorization property 
for the particle hole scattering. This result is consistent with recent calcula- 
tions for the bound state of the massive Thirring model [15-17]. The Bethe 
ansatz as well as the light cone calculations show that there is only one 
bound state on the contrary to the prediction by the S-matrix factorization 
equation. 

Here, we want to discuss the implication of the present result in connection 
with the Yang-Baxter equation. The Yang-Baxter equation is a matrix 
equation which is obtained by imposing certain conditions. The present 
study does not intend to check the Yang-Baxter equation itself. Instead, 
we have only examined the factorization property of the S-matrix at the 
elementary level. Therefore, we have never questioned whether the require- 
ment of the Yang-Baxter equation is physically reasonable or not. We think 
that the requirement of the Yang-Baxter equation is indeed reasonable if 
the S-matrix is factorizable at the elementary level. 

Now, the question is how we can interpret the results (the bound state 
spectrum) predicted by the Yang-Baxter equation when the S-matrix fac- 
torization is violated at the elementary level. In the massive Thirring / sine- 
Gordon model, the situation is now clear. Namely, the spectrum obtained 
by the Yang-Baxter equation with a violation of the S-matrix factorization 
at the elementary level is semiclassical. 

However, this point is only proved for the massive Thirring model. We do 
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not know whether there might be some examples in which the Yang-Baxter 
equation gives exact spectrum even though the S-matrix factorization is 
violated at the elementary level. 
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Figure captions : 



Fig. la : We show the calculated values of -Dipi/i and D2p2h as the function of p^^. 
This is the excited state with — 1, = 1 and rii^ — —1. The solid line 
is drawn with eqs.(5.24) with the parameters in Table 3. 

Fig. lb : The same as Fig. la. This is the excited state with = 2, n^^ = 1 and 
= 2. 

Fig.lc : The same as Fig. la. This is the excited state with = 2, n^^ = 1 and 
rii^ = -2. 
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Table 1 



Rapidities 









at 








101 


0.859 


0.801 


1.142 






401 


0.223 


0.202 


0.247 


Hi = 


5 


1601 


0.0559 


0.0504 


0.0574 






6401 


0.0140 


0.0126 


0.0141 






12801 


0.0070 


0.0063 


0.0070 






101 


1.578 


1.546 


1.863 






401 


0.443 


0.424 


0.486 


Hi = 


10 


1601 


0.1118 


0.1064 


0.1147 






6401 


0.0280 


0.0266 


0.0282 






12801 


0.0140 


0.0133 


0.0140 






101 


2.147 


2.128 


2.390 






401 


0.656 


0.640 


0.713 


Hi = 


15 


1601 


0.1676 


0.1623 


0.1718 






6401 


0.0420 


0.0406 


0.0422 






12801 


0.0210 


0.0203 


0.0210 






101 


2.603 


2.592 


2.806 






401 


0.862 


0.848 


0.926 


Hi = 


20 


1601 


0.2232 


0.2181 


0.2288 






6401 


0.0559 


0.0546 


0.0563 






12801 


0.0280 


0.0273 


0.0281 



We plot the rapidity values for several cases of rii for the vacuum, Ip—lh 
state with rij,, = and 2p — 2h state with rii^ — 1, rii^ — —1. The density 
Po is fixed to po = 20. 
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Table 2 



(^0 = 6.28) 





N 


Dipihinio) 


Eipihiriio) 


D2p2h{'>T'n,ni^) 






101 


0.277 


1.00 xlO"^ 


0.767 


7.31 xlO-'^ 




401 


0.261 


2.40 xlO"^ 


0.764 


1.73 xlO-''^ 


Po = 20 


1601 


0.257 


5.94 xlO"^ 


0.762 


4.25 xlO"^ 




6401 


0.256 


1.48 xlO-^ 


0.763 


1.06 xlO-^ 




12801 


0.256 


7.39 xlQ-' 


0.762 


5.28 xlO-'' 



We plot the values of Dipih, -E^ipih, D2p2h and £'2^2/1 for several cases of 
the number of the particles N for the lowest scattering states { rii^ — 1, 
rii^ — 1, rii^ = —1. ). The density po is fixed to po = 20. 
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Table 3a 



Dipihiuig) {go = 6.28) 





c(i) 


£,(1) 


riio = 1 


-0.359(16) 


0.277(1) 


= -1 


0.359(16) 


-0.277(1) 


^io = 2 


-0.329(14) 


0.275(1) 


riio = -2 


0.329(14) 


-0.275(1) 



We plot the values of the parameters C^-^\ D^^\ The numbers in 
parenthesis indicate error bars in the fitting. 
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Table 3b 



D2p2h{n 


il,'T'i2) (50 = 


= 6.28) 






C(2) 


£,(2) 


1^2 




0.404(8) 


0.626(2) 


22.0 




0.115(1) 


0.454(1) 


12.0 




0.388(5) 


0.603(1) 


18.0 


^ii 1)^12 


-0.404(8) 


-0.626(1) 


22.0 


71 jj^ — 2, 71 — 1 


0.115(1) 


0.454(1) 


12.0 


Tijj^ — 2, 7lj2 — 1 


-0.388(5) 


-0.603(1) 


18.0 



We plot the values of the parameters C*^^\ D^"^^ and K2- The numbers in 
parenthesis indicate error bars in the fitting. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9705140v2 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9705140v2 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9705140v2 



